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The nonadiabatic geometric quantum computation may be achieved using coupled low-capacitance 
Josephson juctions. We show that the nonadiabtic eflFects as well as the adiabatic condition are very 
important for these systems. Moreover, we find that it may be hard to detect the adiabatic Berry's 
phase in this kind of superconducting nanocircuits; but the nonadiabatic phase may be measurable 
with current techniques. Our results may provide useful information for the implementation of 
geometric quantum computation. 
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Quantum computation is now attracting increasing in- 
terest both theoretically and experimentally. So far, a 
number of systems have been proposed as potentially 
viable quantum computer models, including trapped 
ions, cavity quantum electrodynamics, nuclear magnetic 
resonce(NMR), etc In particular, a kind of solid state 
qubits using controllable low-capacitance Josephson juc- 
tions has been paid considerable attention ||,^,^,||. A 
two-qubit gate in many experimental implementations is 
the controlled phase shift, which may be achieved using 
either a conditional dynamic or geometric phase. A re- 
markable feature of the latter lies in that it depends only 
on the geometry of the path executed [||, and therefore 
provides a possibility to perform quantum gate opera- 
tions by an intrinsically fault-tolerant way 1^,^ . 

Recently, several basic ideas of adiabatic geometric 
quantum computation by using NMR S , superconduct- 
ing nanocircuits |^ or trapped ions ^ were proposed. 
However, since some of the quantum gates are quite sen- 
sitive to perturbations of the phase factor of the computa- 
tional basis states, control of the phase factor becomes an 
important issue for both hardware and software. More- 
over, the adiabatic evolution appears to be quite special, 
and thus the nonadiabatic correction on the phase shift 
may need to be considered in some realistic systems as it 
may play a significant role in a whole process ||l^,|ri|,|l2| . 
In this paper, we focus on the nonadiabatic geometric 
phase in superconducting nanocircuits. We indicate that 
the adiabatic Berry's phase, as well as the single qubit 
gate controlled by this phase, may hardly be implemented 
in the present experimental setup. On the other hand, 
since the 2-qubit operations are about 10^ times slower 
than the 1-bit operations [||, the conditional adiabatic 
phase is extremely difficult to be achieved. A serious 
disadvantange of the adiabatic conditional phase shift is 
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that the adiabatic condition requires that the evolution 
time must be much longer than the typical operation time 
To {— fi/Ej with Ej as the Josephson energy), which 
leads to an intrinsical time limitation on the operation 
of quantum gate. Therefore, a generalization to nona- 
diabatic cases is important in controlling the quantum 
gates. We find that the nonadiabatic geometric phase 
shift can also be used to achieve the phase shift in quan- 
tum gates. 




Fig.l 

FIG. 1. Schematic diagram of a quantum computer. The 
j'-th qubit and its probe circuit are displayed in detail. 

We first consider a single qubit using Josephson juc- 
tions described in Ref. (|] ( see the j-th qubit in Fig.l). 
The qubit consists of a superconducting electron box 
formed by an asymmetric SQUID with the Josephson 
coupling El and E2, pierced by a magnetic flux $ and 
subject to an apphed gate voltage 14 = len^jCx, (here 
we omit the subscript j, and 2en^ is the offset charge). In 
the charging regime (where £'1^2 are much smaller than 
the charging energy Er_h) and at low temperatures, the 
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system behaves as an artificial spin- 1/2 particle in a mag- 
netic field, and the effective Hamiltonian reads |l3[ 



(1) 



where crx,y,z are Pauli matrices, and the fictitious field 
B = {Ejcosa, -Ejsina, Echil - 20} (2) 



with Ej = y/{Ei - E2Y + AEiE2Cos'^{'k^/^q), tana = 
{El - E2)tan{-K<^/^o)/{Ei + E2), and $0 = h/2e. In 
this qubit Hamiltonian, charging energy is equivalent to 
the Bz field whereas the Josephson term determines the 
fields in the xy plane. By changing Vx and $ the qubit 
Hamiltonian describes a curve in the parameter space 
{B}. Therefore by adiabatically changing H around a 
circuit in {B}, the eigenstates will accumulate a Berry's 
phase 7b = where the signs ± depend on whether 

the system is in the eigenstate aligned with or against 
the field The solid angle Q, which represents the 
magnetic field trajectory subtends at B = 0, is derived 



BxdtBy — BydfBx 
|B|(i3. + |B|) 



dt. 



(3) 



under the condition B(r) = B(0). 

However, the adiabatic evolution is quite special, and 
thus the generalization to nonadiabatic noncyclic cases 
is of significance. We now recall how to calculate the 
Pancharatnam phase. For a spin-1/2 particle subject 
to an arbitrary magnetic field, each spin state {ip) ~ 
[e-^v/'^cos{0/2), e"^/2sm(6'/2)]^ may be mapped into a 
unit vector n = (sinOcosip, sinOsimp, cos6), with n e a 

unit sphere 5^, via the relation n = {ipl a\ ip), where T 
represents the transposition of matrix. By changing the 
magnetic field, the evolution of spin state is a curve on 
from an initial state (0^, ipi) to a final state (Of^ipf), and 
the Pancharatna m p hase accumulated in this evolution 
was found to be [ [ll| 



function of time t with the period r. We consider the pro- 
cess where a pair of orthogonal states \tp±) evolve cycli- 
cally (but not necessary adiabatically). This process can 
be realized in the present system. Noting that the adi- 
abatic approximation is merely a sufficient but not nec- 
essary condition for the above cyclic evolution, we here 
focus on a nonadiabatic generalization. In this evolution, 
the initial state is given by 



where 



\i^+i0,v)) = [e~'^/^cosi9/2), e'^/^sin{9/2)f, 
\ip-{9,ip)) = [-e-"^/2gm(6l/2), e''^^^cos{0/2)f , 
a+ = cos[{r] — 0i)/2]cosipi/2 — icos[{ri + 9i)/2]sinipi/2, 
a_ = sm[(?/ — 0i)/2]cos(pi/2 + isin[(ri + 0i)/2]sinipi/2, 

with tanr] = Ej{<S? = 0)/Eah, tan0, = [Ej{t)/ Bz{t)]\t^o, 
and tarnpi = [By{t) / Bx{t)]\t=[)- A phase difference be- 
tween I?/; _!_/_) can be introduced by changing H. The 
phases acquired in this way will have both geometrical 
and dynamical components. But the dynamical phase ac- 
cumulated in the whole procedure may be removed , 
thus only the geometric phase remains. By taking into 
account the cyclic condition n(0) = n(r) for l-^i), the 
final state in this case is given by Q 

|^/) = a+e^^|V+(0.,V.))+«-e"'^lV'-(^.,^.)), (5) 

where 7 can be calculated from Eq.(^. The contribution 
from the second term of Eq.(^ vanishes simply because 
n(0) = n(r). Thus the geometric phase considered here 
is the cyclic A A phase. The probability of measuring a 
charge 2e (71 = 1) in the box at the end of this procedure 
is derived as 



Pi — \aj^sin-^ + a^cos-^e ^^t]-^. 
This probability can be simplified to 



(6) 



1 /■ , , sinUpf — ipi) 

7 = -- / [l - cos0)dip + arctan j- -, 

^ Jc cot^cot-j^ -t- cosyipf — ipi) 

(4) 

where C is along the actual evolution curve on S'^, and 
is determined by the equation: dtn{t) — — B(t) x n{t)/h. 
This 7 phase recovers the Aharonov-Anandan (AA) 
phase (Berry's phase) in a cyclic (adiabatic) evolution 

At this stage, we propose how to detect the nonadi- 
abatic or adiabatic geometric phase in the charge qubit 
system. The system is prepared in the ground state of 
the Hamiltonian at ri| = and $ = 0, and then changes 
to the fictitious field B($(t), n^(t)), which is a periodic 



Pi = [1 — cosij] — 0i)cos0i + sin(rj — 0j)szn0jCos27]/2. 

(7) 

when $(0) — 0. Note that Eq.(^ recovers sin'^{j) in 
Ref. even in a nonadiabatic but cyclic evolution [ p^ . 
Thus the nonadiabatic phase may be determined by the 
probability of the charge state in the box at the end of 
this process. It is worth pointing out that the parameters 
T] and 0i in Eq.(^)(or(0)) are fully determined by the 
experimentally controllable parameters and n|, as in 
the adiabatic Berry's phase case 

It is remarkable that the probability obtained in 
Eq.(^(or(|^)) may be directly detected by the dc current 
through the probe junction under a finite bias volt- 
age Vd B. Assume that we have achieved one SQUID 
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qubit as well as the detector circuit, as shown in Fig.l. 
By changing Vjx and $j in time [0, r], the system os- 
cillates between |0) and |1), and the final state would 
be determined by the geometric phase. The measurable 
dc current through the probe junction formulates by the 
processes: |1) emits two electrons to the probe, while |0) 
does nothing. Consequently, the probability described 
by Eq.(||)(or (|^)) as well as the geometric phase may be 
detected by the dc current. 

The single qubit gate may be realized by this geometric 
phase. For example, it is straightforward to check that 
the unitary evolution operator defined by 1-0/) = t/f^lV'i), 
is given by 



cosj isiwy 
isinj cosj 



(8) 



when 9i = and ipi = 0. Clearly, the operation depends 
on the geometric phase 7; 7 = 7r/2 and 7 = 7r/4 procduce 
a spin fiip (NOT-operation) and an equal-weight super- 
position of spin states, respectively. On the other hand, 
the phase-flip gate t/^* = ea;p(— 2i7|l)(l|) (up to an ir- 
relevant over phase ) is derived by 6i — and (pi = 0. 
The noncommutable Ui'' and gates are the two well- 
known universal gates for single-qubit operation. The 
Berry's phase may be used to achieve intrinsical fault- 
tolerant quantum computation since it depends only on 
the evolution path in the parameter space. The nonadia- 
batic cyclic phase is also rather universal in a sense that 
it is the same for a infinite number of possible ways of 
motion along the curves in the projective Hilbert space 
pof . Consequently, the nonadiabatic phase may also be 
used as a tool for some fault-tolerant quantum computa- 
tion. 

We now illustrate how to achieve the cyclic state 
for quantum gates in two processes. The parameters 
(<i>(t), 7i|(i)) in process I change as 
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The path in the parameter space {B} swept out in this 
case is exactly the same as that proposed in Ref. 
Since the evolution in this process is cyclic only under 
the adiabatic condition, we need to answer a key ques- 
tion: whether the adiabatic approximation is valid for 
the given parameters? As for process II , the parameters 
(4>(i),n|(t)) change as 



(10) 



The fictitious field described by Eq.(|l^) guarantees that 
the angle xo = arctan[Ej / {Bz{t) — huj)] (and n^) is time- 
independent. It is found that the state described by the 



vectors n{xo, —cut) in this process evolves cyclically with 
period t = 2Tr/Lu and the A A phase for one cycle is 
given by 7 = 7r(l — cosxo), which may be used to achieve 
the mentioned single-qubit gates geometrically. For the 
present system, the dynamic phase can be removed by 
simply choosing u; = — 4(£'i + E2)Ek[~^EiE2/ {Ei — 
-£'2)^]/7rsm(2xo) with Ek{x) the complete elliptic inte- 
gral of the first kind. 




l=50x„ 
•t = 200i:„ 
i: = 500t. 
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FIG. 2. The trajectories and Bz versus time in process 
I for $,„ = 0.25, n%m = 0.20, E2 = 4Ei = 6.25^lev, and 
Ech = 5.0(^1 4- S2). 

The nonadiabatic effect should be important if r is 
not short. We first consider the evolutions described by 
Eq.(§). Figure 2 shows n^(t) and B^it) = Bz{t)/\B{t)\ 
versus time, with the parameters being the same as 
those in Ref. Q. The deviation of n{t) from B{t){= 
B(t)/|B(t)|) indicates clearly whether or not the adia- 
batic approximation is valid because n{t) almost follows 
the trajectory of the magnetic field B(t) under this ap- 
proximation. It is seen from Fig. 2 that the adiabatic ap- 
proximation is satisfied in the first case when r > SOOtq, 
where tq — h/{Ei + E2) ^ SAps. The adiabatic condi- 
tion for process II is in the same order of magnitude (see 
Fig. 3). It is worth pointing out that the coherence time 
achieved in a single SQUID is merely about 30 ~ 40to 
Q, which is not long enough for the adiabatic evolution, 
implying that the adiabatic condition is not satisfied in 
the above two processes for realistic systems. But, for- 
tunately, the nonadiabatic phase can be measured and 
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used in achieving geometric quantum gates. 

Conditional geometric phase accumulated in one sub- 
system evolution depends on the quantum state of an- 
other sub-system, which may be realized by coupling ca- 
pacitively two asymmetric SQUIDS (see any neighboring 
pair of qubits in Fig.l.). If the couphng capacitance Cy 
is smaller than the others, the Hamiltonian reads 



N 



N-1 



(11) 



where Hi refer to the uncoupled qubits defined in Eq. (|^) 
and H,,,+i = - K,i){nt+i - with 

Ei^i+i = EchCi,i+i/C The gate voltage and mag- 
netic flux can be independently fixed for all qubits. We 
address firstly a two-qubit operation, e.g., i and j qubits 
are two neighbour qubits with the i-th as the control 
qubit and the J-th as the target qubit. The fictitious field 
on the target qubit is [Ej{<^j)cosaj,—Ej{<^j)sinaj,B''^] 
with i?^ = Ech{^ — 2n'^ j)+Eij{n'^ i—l), where I represent 
the control qubit state or I. Obviously, the geometric 
phase 7j for j-th qubit in decoupled case is different from 
7^ even the changeings of n| j) are the same, where 
7j is the geometric phase of the target qubit when the 
charge state of the control qubit is I. 7] may be di- 
rectly derived from Eq.(^. It is worth to pointing out 
that the state described by the vector n{x'',—^t) with 
Xg = arctan[Ej/{B'^ — hu)] is still a cyclic evolution, 
and may be used to achieve the two-qubit operation. In 
terms of the basis {|00), |01), |10), the unitary op- 

erator to describe the two-qubit gate is given by 



U, 



(7°,7j) 



diag{e , e^'^'j 



(12) 



The combination with single-bit operations allow us 
to perform the XOR gate. The unitary opera- 
tion for XOR gate can be obtained by Uxor — 
[I ® U''^{tt/A)] t/(o,;^) [/ «) J7i'^(7r/4)]^ with / as a 2 x 2 
unit matrix. This XOR gate together with single qubit 
gates constitutes a universality: they are sufficient for all 
manipulations required for quantum computation [|l8| . 
Therefore, all the elements of quantum computation may 
be achievable by (nonadiabatic) geometric phase. 

We now compute the geometric phases required by 
the spin flip operation (7 — 3tt/2) and NOT operation 
(7 — 37r/4) accumulated in the second process. The com- 
parison of the nonadiabatic geometric phase 7 with is 
shown in Fig. 3, where the 7a is the phase calculated un- 
der the adiabatic approximation. It is seen that the 7^ 
deviates evidently from the 7 for r < ISOtq. Thus the op- 
eration time required by the adiabatic condition in both 
processes I and II is in the same order of magnitude. 
The dynamic phases can be removed when t ^ 3.57to 
for 7 = 37r/2 and t ^ 4.I2to for 7 — 3tt/A, respecitively. 
Therefore, by accurately controlling the parameters 
and n^, we may control the state in the projective Hilbert 



space. It is striking that the present operation time 
for nonadiabatic geometric gates is much shorter than 
that with the adiabatic scheme. Note that the coherence 
time achieved in a single SQUID by current technology 
is about 30 ~ 40to , implying that tens of geometric 
NOT operation may be achieved experimentally. There- 
fore, the generalization of the adiabatic phase to the 
nonadiabatic case is of significance since the coherence 
time achieved in charge qubit in Josephson Junctions is 
short. Moreover, the large number qubits required for 
useful computation may be devised by a network similar 
to Fig.l. 
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FIG. 3. The comparison of adiabatic phases with nonadia- 
batic ones. The parameters are the same as those in Fig. 2. 

In conclusion, we study how to detect the nonadiabatic 
phase in superconducting nanocircuits, and the possibil- 
ity to use the nonadiabatic phase as a tool to achieve the 
quantum computation. 
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